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I. INTRODUCTION 



Chiral Particle Decay of Heavy-Light Mesons 
in a Relativistic Potential Model 

Takayuki MatsukiQ and Koichi Seo^j] 

[ Tokyo Kasei University, 1-18-1 Kaga, Itabashi, Tokyo 173-8602, Japan 
2 Gifu City Women's College, 7-1 Hito-ichiba Kitamachi, Gifn 501-0192, Japan 
(Dated: December 30, 2011) 

Partial decay widths of the heavy-light mesons, D,D S ,B, and B s , emitting one chiral particle {it or K) are 
evaluated in the framework of a relativistic potential model. Decay amplitudes are calculated by keeping the 
Lorentz invariance as far as possible and use has been made of the Lorentz-boosted relativistic wave functions 
of the heavy-light mesons. One of predictions of our calculation is very narrow widths of a few keV for yet 
undsicovered B,(0 + , 1 + ) mesons corresponding to 2S + l Lj = 3 Po and " 3 Pi" assuming their masses to be 5617 
and 5682 MeV, respectively, as calculated in our former paper. In the course of our calculation, new sum rules 
f"^ ' are discovered on the decay widths in the limit of mg — > oo. Among these rules, r(£)* (2317) — > D s + 7t) = 

(N ■ T(Z) jl (2460) -> D' s + it) and r(s^ (5615) -> B s +7t) = r(S sl (5679) -> S* + it) are predicted to hold with a 

' good accuracy. 

<D ' 

Q ■ PACS numbers: 13.25.Gv, 14.40.Pq, 13.75.Lb 

(N 

We have been trying to explain the mass spectrum of the heavy-light mesons, including the famous D s o(2317) and 
^"'i D*j(2460), by a relativistic potential model with a linear potential and the Coulombic potential. In our model besides the 
<-j | current quark mass m q (m Ut d ~ 10 MeV) the constant term b is introduced in the scalar potential. They are independent parame- 
'— 1 , ters in general. Howevere, if we restrict our computation up to the first order in 1 /mg, the current quark mass m q and the constant 
b appear only in the combination m q + b. Hence we set b — and call m q in place of m q + b. In this model, we have successfully 
reproduced many of the experimental mass spectrum of D, D s , B, and B s with a fairy good accuracy including radially excited 
states!!. 

If} [ To confirm the validity of our relativistic potential model, we have calculated the semileptonic weak form factors (Isgur-Wise 
functions) for the process B — > D^€v in Ref. J2l to obtain reasonable results compared with the experiments. In the present 
paper [130 , we calculate the decay processes of the heavy-light mesons with one chiral particle emission. Although the decay 
rates of these processes have been examined by Goity and Roberts [0] and by Di Pierro and Eichten J^t] in a relativistic potential 
model, we would like to see the effects of using our relativistic wave functions as well as Lorentz covariant amplitudes that have 
been neglected in their papers. The authors of Refs. Q] and Jit] have used wave functions of the heavy-light mesons in the rest 
frame to calculate transition amplitudes. That is, they have neglected not only the relativistic effects of the wave functions but 
also the recoil effects of the heavy-light mesons. We claim that this recoil effect is essential to calculate the relativistic transition 
amplitude because the plane wave function of the emitted particle is replaced by another phase factor, which involves the heavy 
quark mass and the velocity of the heavy-light meson. The plane wave, which, for instance, is used in a hydrogen atom, is 
inserted without any criticism, which may lead to the erroneous conclusion. In the course of our calculation, a couple of new 
sum rules on the decay widths are discovered, which hold among the different decay processes in the heavy quark limit, i.e., 
niQ — > oo. The breaking of the sum rules is of the order of (M, - M/)/mg, which is at most 10 % or so with the initial and final 
heavy-light meson masses, M; and Mj. 

In Sec. |nj the wave function of the heavy-light meson in the moving frame is related to the one in the rest frame by regarding 
the heavy quark to be at rest in the heavy-light meson. Then we derive the basic formula to express the transition amplitude in 
terms of the wave functions of the heavy-light meson in the rest frame. The plane wave of the emitted particle moving in the 
z-direction, exp(-jfe), used in other works 0] and JH] is found to be replaced by another phase factor, exp(-2/mgVz), with V 
being the heavy-light meson velocity in the Breit frame. In the limit of ihq — > oo, the heavy quark mass ihq coincides with the 
heavy-light meson mass and these two phase factors become equal to each other. 

Interaction between light quarks and chiral particles {it or K) is assumed a la Georgi-Manohar |6]. The matrix elements of 
axial-vector currents are dissolved into several form factors, which are expressed in terms of the radial wave functions. 
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In Sec. |nil we numerically evaluate the matrix elements of axial-vector currents. Except for the axial- vector coupling constant, 
We adopt the parameters of the model which have been determined in the previous paper to fit with the experimental mass spectra. 
The axial-vector coupling constant is adjusted so that the theoretical decay widths of D* + (2010) coincide with the experimental 
values, whose derivation is still contentious and is discussed in Sec. [V] We predict the decay widths of the various decay modes, 
most of which have not yet been experimentally observed. The total widths of some heavy-light mesons are experimentally 
known, which are compared with our theoretical results. Because D s (1~), D s (0 + ), and D s (1 + ) are not allowed kinematically to 
emit one K meson, assuming the mixing of n° and 77 we have the narrow decay widths of these particles. Likewise, we predict 
very narrow decay widths of B s (0 + ) and B s (1 + ), which are not yet observed and hence we have used our theoretical values of 
their masses, 5615 MeV and 5679 MeV, respectively. 

In Sec. |IV1 new sum rules on the decay widths are presented in the limit of m Q — > 00. In this limit the spin of the heavy quark 
is decoupled from the orbital angular momentum and the spin of the light quark, and the heavy quark is no more than a spectator 
in the decay process. The radial wave functions of J p and (J+ l) p coinside with each other and their masses become degenerate. 
We call that these degenerate particles belong to a spin multiplet. The decay widths of particles in the same multiplet become 
equal if we sum over the final states in another multiplet. 

Sec. [V] is devoted to the conclusions and discussion. In Appendix A, the matrix elements of the axial-vector currents are 
described by the polarization vectors/tensors in the rest frame of each heavy-light meson. In Appendix B, the properties of 
the polarization vectors are given as well as the explicit forms of angular-spin part of the wave function. In Appendix C, the 
transition amplitudes are expressed in terms of the radial wave functions and the spherical Bessel function, f(x) = sinx/x. 



II. LORENTZ INVARIANT EVALUATION OF TRANSITION AMPLITUDES 



The wave function of a heavy-light meson with a finite momentum is defined as 



<O\qiat)Q0,t)\P> = < 1 qfitf- X ? , 0) Q0- A f , 0) I P > e" 



(1) 



where X% (= gx + (1 - denotes the position of the heavy-light meson and £ is a free parameter. If we set £ = or 1, then X^ 
coincides with the position of the heavy quark or light quark, respectively. 

We assume chiral interaction of pseudoscalar mesons {n and K) with light quarks. In the present paper we only compute one 
pseudoscalar particle emission from heavy-light mesons, and then the relevant interaction Lagrangian is as follows: 



■Cint — 



V5/„ 



(2) 



where g is a dimensionless coupling constant and f„ is the pion decay constant. Here the flavor SU(3) symmetry is assumed and 
<f>i/s represent the octet meson fields, that is, 



M = ^2 



( Ti^ + 

71 

K 



K + 
K° 



(3) 



The mixing of 7r° and 77 is taken into account with a small parameter e as follows: 



VT 



+ e- 



1 e 
-e 1 



(4) 



The transition amplitude for a heavy-light meson with a momentum P into another heavy-light meson with P' and one 
pseudoscalar particle with k is transformed by inserting the number operator of the heavy quark and by neglecting the sea 
quark effects as follows: 



Mr 



= J d 4 x <P',k\£ int (x)\P> 

= jd 4 x < P', k I £ int {x) J d 3 y Q\(y)Q e {y) \ P > 

* J d 4 x J d'y { j < P' \Q](y)qf (J01 0x0 \q c j(x)Q e (y)\ P > e ik x 
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(2n) 4 6 4 (P -P'-k) J d 3 z tr P') O u P)\ e~ iU (5) 

(2*) V(P - P' - Jfc) J ^ztr^W^POOy^^/O] 



(6) 

In Ref. J2l the following approximate relation has been obtained between the wave function with a finite momentum (P3 = 
yMV) and the one at rest, 

0g(jf, P) * GoyG^Cxx, yz; M) e ^-m e )v rZ) (7) 

and the following relation is derived likewise, 

^(f; P) * GayGfxi/fysixx, yz; M) e -'" n ^\ (8) 

where G is the boost operator from the rest frame to the moving frame and ty(x; M) is the wave function in the rest frame. The 
Lorentz transformation from the rest frame to the moving frame produces the time difference between two constituent quarks, 
and this time difference is compensated by the free propagation of the heavy quark, which gives rise to the phase factors in 
Eqs.© and ®. 

If we evaluate the transition amplitude in the Breit frame, where the initial heavy-light meson is moving to the z direction with 
a velocity V and the final heavy-light meson is moving in the opposite direction with the same velocity V, and insert Eqs.((7]i and 
© into Eqs.© and ©, we obtain the following unique expression, 

AW< ~ (2n) 4 6 4 (P - P' - k) 

x j d\ tr [(G" 1 ) VJOtj., yz; M')G- { O u G ifr je (x ± , yz; M)G T ] e ~ 2im a? Vz 

= (2tt) V(P -P'-k) 

xy- 1 j d\ tr [(G- l ) T ^ e {x; M^G' 1 O u G if/ je (x; M)G T ] e~ 2im a Vz . (9) 

Here two points should be noted. First this is independent of £ as it should be. Next the the phase factor e~ 2 " neVz has appeared 
in place of the wave function of the emitted pseudoscalar particle, e~' kz , that is absorbed into the delta function to preserve the 
four-momentum conservation in either Eq.© or ©■ This is the result of taking into account of the heavy-meson's recoil effect. 

In general the matrix elements of the axial-vector current between the various spin states of the heavy-light mesons have the 
following tensor structures, 

= (1 + (o)e\^ A Uw) 

+(si ■ v 2 ) ((v, + va)^(a») + (v, - v 2 )^(w)) (k = 1, 10), (10) 



< 0-|^|0 + > 



< i-\h„\o + > 

VM2M1 

<o-\j 5 ,\i + > 



(VI + V2)^l(co) + (vi - v 2 )^(oS), (11) 
e^vy^^ico), (12) 



VM2M1 

< n/sjr > 

i VM2M1 

< 0-|75^[2 + > 
/ \/M 2 M\ 

< Hi5 M i2 + > 



= e u 



r v^ef^f(w) (* = 4,5), (13) 



(£2 • si)(vi + v 2 )^*{(cj) + (e 2 • ei)(vi - v 2 ) M ^(w) 

+(4 ■ vi)£i„O w ) + ( £ i ■ V2)fi^(w) (* = 6, 7), (14) 
ej^v^Cw) + e^v^ (v]^g(w) + v 2/J ^ 8 3 V)} , (15) 

{(V! + viYC^) + (vi " V^C^)) 



VM 2 Mi 
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< 11^11- > 
i ^M 2 M\ 



c /jvpo~ 



Wssfv^sfvlvl {(vi + v^^V) + (vj - v 2 )^^M} , 
e!4 P ((vi + virffiW + (vi - VaT&V)} 

+e of} y S e a l efv y 1 V2 {(vi + v^^^ici) + {v\ - v 2 )^\co)] , 



(16) 



(17) 



where the initial and final quantities have subindices 1 and 2, respectively, and s- lllv are polarization vector and tensor, vi and 
v 2 are velocity vectors defined by P\ = M\\'\ and P 2 = M 2 v 2 , which in the Breit frame have components, v° = v° = y, v\ = 
— vf = •yV, and <w = vi • V2 = Pi ■ P 2 I{M\M 2 ), which in the Breit frame is related to V by 

1 + V 2 



co = r 2 (l + V 2 ) 



l-V 2 ' 



(18) 



In our model the 1 + mass-eigenstates appearing in Eqs. (13[ and ( TBI are characterized by the quantum number of j = L + s q , and 
3 P\ and l Pi dominant states are denoted by " 3 Pi" and " l P\", respectively. This is the reason why we have two kinds of form 
factors in these equations, i.e., k — 4, 5 and k = 6,7. There are also two 1" states with k = 1, 10 in Eq. ( fTOb corresponding to 
3 S i and 3 Di states, respectively. Each componet of the matrix elements of the axial-vector current in the Breit frame is written 
down with the polarization vectors/tensors in each rest frame of the initial and final mesons in Appendix A. 

In order to obtain explicit forms of the form factors £a's, we need to calculate the matrix elements by inserting the explicit 
forms of the wave functions. According to Eq.® the matrix elements of the axial vector current are evaluated in terms of the 
wave functions in the rest frame of the heavy-light mesons, fci?, M), as follows: 



2y/MiM 2 

«M?;M) 



y 



\ d 3 x — Ur-tr 

J 4nr 2 2 



y* 2 (m 2 , i(n ■ cf)v 2 ) G l (p u o- u o- 2 ,o-T)G 



Ml 

—i(H ■ <t)v\ 



iyi 



ia(r) 
-On ■ cf)vi(r) 



yu 



(19) 
(20) 



where q = -2mgV, <x,- = cr, ® l 2x2 , and ft = r/r. By substitution of the explicit expression for G, the matrix elements of each 
component of the axial vector current are obtained as follows: 



2 y]M\M 2 

< ^2lj5 3 |^l > 



2 VM,M 2 



2 VMiM 2 



/l 1 
j3x ^2"( M2Vl - v 2 wi)-tr[y|(« • (?)yi] e' qz , 

/I i 
d 3 x-^-jU 2 ui-ti\y* 2 o- 3 yi] e lqz 

/I i 
d 3 x ^— j v 2 vi -tr [y 2 (« • oV 3 (« • <?);yi] e ,V/ \ 

foil 

I d x-^—^u 2 ui-tT[y* 2 o-iyi] e ,qz 

foil 

+ I c/ x— ^v 2 vi-tr[y;(n • ff)cr0 ■ &)yi] e' q ~ 

/l | [ 

d 3 x ^2 e ^i V 2 tX y * 2 { v 2"i(« ' (? )°"j - u 2 vio-j{n ■ (?)}yi 



(21) 



(22) 



(/= 1,2). 



(23) 



The angular-spin part of wave functions, y,(f2), are given in Appendix B. 

The form factors £a's are expressed in terms of the radial wave functions u's and v's in Appendix C, where we have adopted 
only the terms up to the first order in V. The formula for the decay width, 



r = 



* l . — y i 



<P 2 ,k\£ inl (0)\P l >\ 2 , 



(24) 



is also evaluated in terms of ^'s in Appendix C, where k^ is the momentum of the chiral particle in the rest frame of the parent 
heavy-light meson. 

If the mass difference of the parent heavy-light meson with strangeness and daughter non-strange meson is smaller than the 
A"-on mass, then the parent heavy-light meson with strangeness decays into a heavy-light meson with strangeness and a neutral 
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pic 
17] 



ion through the n Q -r\ mixing. The mixing parameter e has been related to the current quark masses due to Gasser and Leutwyler 



as follows: 



^ md - m " -l.OOxlO- 2 (25) 



4 m s - (m u + md)/2 



As a result, the formula of the decay width is multiplied by the factor, e 2 (= 1 .00 x 10 4 ), in addition to the factor 2/3 (the square 
of the coefficient of </>g in the octet meson matrix (01). 



III. NUMERICAL RESULTS 



In the leading order of l/mg expansion, the radial wave functions u(r) and v(r) are determined by the following Dirac 
equation, 

± + h _ mq _ S +V)\v(r)) h \v(r)}> {Zb) 

where S and V are confining scalar and Coulombic potentials parametrized by S (r) — b + and V(r) — ^f. k is the eigenvalue 

of the operator -B q {L ■ & q + 1). The states with the total angular momentum j and j + 1 have the same value of k, and they have 
the same radial wave function and the same eigenvalue E. For example, S -wave states 0~ and 1~ with k — -1 are degenerate. 
Taking into account the asymptotic behaviors at infinity and at the origin, u(r) and \>(r) are approximated by the following forms, 

7 

u{r) = r 7 exp{-r 2 /a 2 -(b + m q )r}Ya ( f ) r i , (27) 

i=0 

7 

v(r) = r y exp{-r 2 /a 2 -(fe + m 9 )r}^]fl|"V, (28) 



where y = ^ 2 - (^-) 2 . The values of the parameters used in the numerical calculation are listed in Table U which are 
determined so as to reproduce the mass spectra of the heavy-light mesons in the order of l/mg. (See Ref. |Q]].) In the present 
calculation the constant b and the light quark mass appear only in the linear combination, and the values of this combination are 
given in the table. 



TABLE I: Values of the parameters 





a 


b + m„(= m d ) 


b + m s 


m c 


m b 


8 2 




(GeV- 1 ) 


(MeV/c 2 ) 


(MeV/c 2 ) 


(MeV/c 2 ) 


(MeV/c 2 ) 




0.259(D) 
0.392(B) 


1.937 


86 


168 


1023 


4634 


0.608 



Here different values of the strong coupling a s are used for D/D s and B/B s mesons. This parameter set is not a unique 
solution, but the results are not so changed even if we take another set of values given in Ref. yfl. The axial-vector coupling 
constant g is determined so that the calculated decay widths of D* + (2010) agree with the experimental values. The value of 
g 2 is taken as 0.608 with a 27% uncertainty. The errors of prediction are caused by the uncertainty of our parameters besides 
g 2 , and only the central values of the calculated widths are shown in Table 2, 3, 4 and 5. k^ is the momentum of the emitted 
chiral particle in the rest frame of the parent heavy-light meson. The observed widths in our tables are taken from Ref. 10], and 
experimental values depending on model assumptions are omitted. 

TABLE II: Numerical evaluation of the decay widths of excited D mesons 



Initial State 


Final 






^exp 




State 


(MeV/c) 


(MeV) 


(MeV) 


D*° 








not allowed 


( 3 Si) 


nO 

D n 


43.1 


0.042 


< 1.3 


jyt± 


D 71 


39.5 


Input 


r\ f\£L C i A A1 o 

0.065 ± 0.018 


( 3 5l) 


D*n° 


38.2 


Input 


0.029 ± 0.08 






414.2 


0.99 x 10 2 




D*,(2400)° 


D°7I° 


419.4 


0.50 x 10 2 






Dn(sum) 
all 




2.5 x 10 2 


261 ±50 




DV 


461.3 


1.2 x 10 2 




05(2400)* 




458.5 


0.61 x 10 2 




( Po) 


D7r(sum) 
all 




1 o * * 1 rv2 

1.8 x 10 


283 + 24+ 34 


D n (2400) 


D n 




u 


not seen 


£>i(2430) 


Dn 













358.7 


0.76 x 10 2 




Z) 1 (2430)° 


D*°7T° 


363.1 


0.39 x 10 2 




(" 3 i°i") 


D*n(sum) 
all 




1.9 x 10 2 


384 +10? ±74 

— /J 


£>i(2420) 


Dn 







not seen 






354.5 


2.7 




D,(2420)° 


D*°7T° 


358.9 


1.5 




C'^i") 


D*7r(sum) 
all 




6.9 


20.4 ± 1.7 






358.4 


2.9 




Di(2420) 


D TT 


357.0 


1.4 






D*7r(sum) 
all 




4.3 


25 ± 6 




D±7T T 


505.5 


6.8 




D*(2460)° 


D°n° 


510.2 


3.6 




(^2) 




389.2 


2.6 






D*°7T° 


393.4 


1.4 






Dn + D*7r(sum) 




24 






all 






43 ±4 




DV 


508.4 


7.1 




Z)*(2460) ± 


D ± 7T° 


505.6 


3.4 




( 3 ^2) 


D*V 


391.2 


2.7 






D^n 


389.7 


1.3 






Dn + D*7r(sum) 




15 






all 






37 ±6 




DV 


739.6 


4.8 




£>(2760)° 


D°n° 


743.6 


2.4 




( 3 £>i) 




638.9 


0.79 






D*V 


642.1 


0.41 






Dn + D*7r(sum) 




14 






all 






60.9 ±5.1 




DV 


742.9 


4.9 




£((2760)* 


D ± 7T° 


740.3 


2.4 




( 3 £>i) 


£rV 


641.3 


0.81 








639.6 


0.39 






Dn + D*7r(sum) 




8.5 
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Initial State 
{ 2S+l Lj) 


Final 
State 


k R 
(MeV/c) 


(MeV) 


r 

A exp 

(MeV) 




all 






no exp 



The narrow width of D*{\) is due to the smallness of the Q value. The width of £>i(2420) is smaller than that of £>i(2430) by 
an order of magnitude. The values of k of the former and the latter are assumed to be -2 and 1, respectively, and their angular 
momenta neglecting the heavy quark spin denoted as j( are | and ^respectively. By the parity and the angular momentum 
conservation, 1 + is able to decay into 1~ and one pseudoscalar particle with an orbital angular momentum of L = or 2. Since 
jc of the daughter meson D*(l~) is i D\{2A20){k = -2) emits a pseudoscalar particle with L — 2 only, while Di(2430)(&: =1) 
emits one with L = 0. This is the widely accepted explanation of the difference of the decay widths of D 1 (2420) and D 1 (2430). 
As for the time componet of the axial vector current, this argument is true and the ratio of the matrix elements is almost one- 
thousandth. As a matter of fact, the space component parallel to the momentum of the emitted particle in the Breit frame is not 
so suppresed, and the ratio of the widths amounts to 10. 

TABLE III: Numerical evaluation of the decay widths of excited D s mesons 



Initial State 


Final 


k R 


F f /, 


r 

A exp 


( 2s+l Lj) 


State 


(MeV/c) 


(MeV) 


(MeV) 


ores,) 




47.8 


8.0 x 10~ 6 


< 0.11 


£>: (2317) ± 


n±^> 


297.7 


3.8 x 10~ 3 


< 3.8 


( 3 Po) 


D^TT 


148.0 







D. sl (2460) ± 


DfTfi 


424.8 







C 3 Pi") 




297.3 


3.9 x 10 3 


< 1.7 




D^K* 


390.4 





not seen 


0,1(2536)* 


D ± K° 


385.6 





not seen 


c i pn 


D *o K ± 


165.3 


0.066 






D »± K o 


159.4 


0.055 






D*K(sum) 




0.12 






all 






<2.3 






433.9 


3.4 




D j2 (2573) ± 


D ± K° 


429.4 


3.3 




( 3 P 2 ) 


D *o K ± 


242.6 


0.27 






D *± K o 


238.3 


0.24 






DK + D*K(s\xm) 




7.2 






all 






20 ±5 




D^K* 


673.2 


6.6 




0,(2818)* 


D ± K° 


669.8 


6.4 




( 3 £»i) 


D *o K ± 


547.3 


1.1 




not observed 


D »± K o 


544.8 


1.0 






DK + D*K(sum) 




15 





D*, Z)* (2317), and D s i(2460) are kinematically forbidden to emit a strange particle, and their decays proceed through the 
7T° - T] mixing. The mixing parameter e is assumed to be 0.01, and these decay widths are smaller than those of the heavy mesons 
allowed to decay with an emission of a strange particle by four orders of magnitude. 

TABLE IV: Numerical evaluation of the decay widths of excited B mesons 



Initial State 


Final 


k R 




r 

A exp 




State 


(MeV/c) 


(MeV) 


(MeV) 


[4] B*(5590)° 


B ± n T 


269.9 


35 




( 3 A>) 


B n° 


271.8 


17 





8 



illlLldl oLaLC 


Final 




J- th 


r 

A exp 


(2S + 1 j \ 
\ Lj ) 


OLaLe 


^ ivie v / c ) 


^ivie v ) 


^ivie v j 


not observed 


Z?7r(sum) 




1 


no exp 


#i(5646) 


r>*± T 

B 7t 


280.7 


37 




C 3p n 


ZTV 


282.8 


19 




not observed 


B*^-(sum) 




93 


no exp 


Bi(5721)° 


B^TT* 


360.0 


7.5 




(" l Pi") 


r>*0 
B U 7T U 

B*7r(sum) 
all 


361.7 


3.8 

1 C\ 

19 


no exp 




B ± n :f 


424.4 


6.4 




Kj(5747)° 


B°n° 


425.5 


3.2 




(^2) 


B^JT* 


379.5 


5.7 






B*°7T° 


381.1 


2.9 






Bn + B*7r(sum) 




30 






all 






227+3.8 +3.2 
-3.2 -10.2 




B ± ii* 


651.0 


5.8 




5(5985)° 


B°7T° 


651.6 


2.9 




( 3 £>i) 


B* ± 7T T 


609.4 


2.3 




not observed 


B*°n° 


610.3 


1.1 






Bn + ZT7r(sum) 




20 


no exp 



3 _Po " 3 -Pi" (1 + ) and 3 Di (1 + ) states of the bottomed meson are not yet established, and their masses are taken from our 
paper {XO - The only observed decay width of B* 2 (51A1) is consistent with our prediction. 



TABLE V: Numerical evaluation of the decay widths of excited B s mesons 



Initial State 


Final 


k R 




I~ exp 


(2S + 1 j \ 

{ Ljj 


State 


(MeV/c) 


(Mev) 


(Mev) 


B s0 (5615) u 


dO 


204.2 


1 £ w 1 A — 3 

l.o x 1U 


no exp 


( Po) 


D*0 









not observed 










B. s i(5679) u 


r>0 









(" 3 fi") 


D*0 

B,r 


224.0 


1.9 x 10~ 3 


no exp 


not observed 














229.2 







B. s i(5830)° 


B°K° 


230.5 









B* ± K T 


93.6 


0.011 






B* u a u 


98.3 


0.014 






B a (sum) 




0.036 






all 






no exp 




B ± K T 


250.6 


0.56 




B; 2 (5840)° 


B°K° 


251.8 


0.58 








135.3 


0.041 






B"K" 


138.6 


0.046 






DA + D A(,sumj 




1 8 
l.o 






all 






no exp 




B ± K T 


523.1 


10 




B s (6025)° 


B°K« 


523.5 


10 





10 



Initial State 


Final 


k R 


T, h 


r exp 


t2S + l t \ 

[ Ljj 


State 


(MeV/c) 


(MeV) 


(MeV) 


( 3 £>i) 




465.9 


3.5 




not observed 


B*°K° 


466.8 


3.5 






BK + B*K(sum) 




41 


no exp 



3 Pq(0 + ), " 3 f 1 "(l + ) and 3 Di(l + ) states of the B s meson are not yet found, and their masses are also taken from our paper []]]. 
As is the D s meson, the narrow widths of 3 Po(0 + ) and " 3 fY'(l + ) states of the B s meson are explained by kinematics. The 
predicted mass values of these particles are below the the threshold of Konic decay, and they decay only through the n° - 77 
mixing. 



IV. NEW SUM RULES ON THE DECAY WIDTHS IN THE LIMIT OF m Q -> co 



In the limit of nig — * 00, the spin of the heavy quark is decoupled from the orbital angular momentum and the spin of the 
light quark, and the radial wave functioin is solved for each eigenvalue of k. The states with J p and (J + l) p take the same value 
of k, and they belong to the same spin multiplet. They are mixed by the spin rotation of the heavy quark, %. If we denote the 
transformation of the state as 



H\P,k,P >=\P,k,q >, 



then 



^\P,k,p >< P, k, p I ft f = Y I P, k, q >< P, k, q \ 



(29) 



(30) 



where P stands for the four-mometum of the heavy-light meson, k is the eigenvalue to discriminate the spin multiplet, and p 
and q label the polarization of the total spin. The sum is taken over all states belonging to the same spin multiplet. Because the 
axial-vector current for the light quark and hence the interaction given by Eq. (f2]are invariant under %, we obtain the following 
relation among the transition amplitudes 



Y^\< P\k',p'\Lin,\PXp >| 2 

p' 

= Yj\< P , ,kf,p'\R f 9i£ba^9<\P,k,p > 



< P,k,pW£in,K 



Y\P',V,p' >< P',k',p' 

< p' 

= < P,k,q\£ int Y J \P',k',q' X P',k',q'\£ int \P,k,q > 
1' 

= Yj\< P',k',q'\£ im \P,k,q >\ 2 - 



1?£ int nP,Kp > 



(31) 



Neglecting the mass difference of the mesons in the same multiplet, we obtain the following relation on the deccay widths, 

Y j r(\P,k,p)^\P',k',p') + 7t/K) = Y j r(\P,k,q)^\P',k',p') + 7t/K), (32) 
p' p' 

where the sum is taken over all states belonging to the same spin multiplet and p and q are in the same spin multiplet also. 



j^-j J] T(\P,k,p{J))^\P',k',p') + nlK) 
J pdip' 

—J— J] r(\P,k,p(j+D)^\P',k',p') + n/K). 
KJ ' p(j+i),p' 



(33) 



We are able to explicitly confirm these sum rules in the case of the initial states (0 + , 1 + ) with k — 1 and (1 + , 2 + ) with k = -2 
and the final states (0~, 1~) with k — -1 by using the expressions of the decay widths given in Appendix C. 
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As a matter of fact, there exist mass difference between the mesons belonging to the same spin multiplet, and the sum rules are 
violated more or less by the kinematical phase factor. Among the various decay modes the following sum rules are not largely 
affected by the kinematical phase factor because of relatively small mass difference between initial and final heavy-light mesons, 
and these rules should be confirmed by experiments, 

T (z)* (0 + , 2317) -> D s (0~, 1968) + n) = T (D sl (l + , 2460) -> D* s (r, 2112) + n) , 

r (fi; o (0 + , 5615) -> B s (0~, 5366) + n) = T (B sl (l + , 5679) -> 5415) + n) . (34) 

Here the decay modes + — » 1~ + n and 1 + — » 0" + n are forbidden and the sum rules are simply reduced to the above equations. 



V. CONCLUSIONS AND DISCUSSION 



The wave function of the heavy-light meson in the moving frame is approximately related to the one in the rest frame. 
Then the the transition amplitudes of the excited states of the heavy-light meson to the lower states by emitting a chiral particle 
are expressed in terms of the wave functions of the heavy-light meson in the rest frame. The plane wave of the emitted particle, 
e~' kz , has been inserted without an explanation by all authors in the preceding works, but we have found that a phase factor 
e -2im g Vz should b e inserted instead of e~' kz . The boost operator affects the matrix elements of the components of the the axial- 
vector currents perpendicular to the momentum of the emitted chiral particle. However these componets are irrelevant in our 
calculation because the matrix elements of the axial-vector current are contarcted with the momentum of the emitted chiral 
particle. 

We have used the wave functions of the heavy-light meson obtained by the relativistic potential model in Ref. [1]. Beacause 
3 Po (0 + ), " 3 Pi" (1 + ) and 3 £>i (1 + ) states of the B and B s mesons have not yet been observed, we have used our predicted mass 
values for these particles [Q]]. 

The partial decay widths of the excited heavy-light meson emitting one chiral particle are numerically evaluated, and the 
predicted values are consistent with the experimentally observed total decay widths of excited heavy-light mesons. 

One of notable predictions is the narrow decay widths of undiscovered B S (Q + , 1 + ) mesons. Our predicted mass values of these 
are below the threshold of the Konic decay, and they decay through the 7r° - rj mixing with narrow decay widths of a few keV. 
Our predictions are awaited for the experimental confirmation. 

The parameters used in our calculation are not much different from those used by Di Pierro and Eichten @1, but the plain 
wave of the emitted particle(e~' fc ) in the transition amplitudes is replaced by another phase factor e - 2,m Q Vz due to the recoil effect. 
This replacement makes the decay widths large in general, because the oscillatory cancellation due to this phase factor becomes 
weak. 

There are numerous studies to compute the decay widths using the Schwinger-Dyson amplitudes, for example, Ref. Js}] and 
references therein. The coupling constant g appearing in these articles and our axial- vector coupling constant g appearing in Eq. 
(fJJ are related to each other in the case of D* — » D + n or kinematically forbidden case of B* — » B + n as follows: 

^ (O'Mjs^') = (Htt\H*) (in their notation), 

left hand side * 2 l^pE± (^U) _ ^(0) ^ 

fn 

_ ig^mm 3 (i, (i, _ m 

f k r i A2 e i' yn M - u), 

Jit 

■ u*u a -a 2gyfWM ii 
right hand side = Ar 63, 

fit 

?7 (1) 's are defined in Eqs. (CI) and (C2) in Appendix C, which are the ovelapping integrals of the initial and final wave functions. 
Coupling constants thus determined may be used in the study of dissociation processes 7r + //(/'—>£) + £) or 7r + Y— > B + Bby 
exchanging D* or B*. We have obtained grf^ = V0.608 x 0.74 = 0.577 while they gave the value g — 0.53 using the leptonic 
decay of the heavy-light system. ^ Agreement of these values is not surprising because the authors of Ref. |@] as well as we 
have fitted couplings with the experimental data. 

In the course of the calculation, new sum rules are found to hold in the limit of ihq — > 00. Among these rules, F(D* (2317) — > 
D s + n) = r(D s i(2460) -> D* + n) and r(B* (5615) -> B s + n) = T(B J i(5679) -> B* + n) are supposed to hold with a good 
accuracy and should be verified by future experiments. 

The radiative decays of the heavy-light mesons are under study in the same formalism as the present one and the results will 
be published in near future. 
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Appendix A: Tensor structures of the matrix elements of the axial-vector current 

The matrix elements of the axial-vector current are expressed by the polarization vectors/tensors in the rest frame of each 
heavy-light meson as follows: 



<oifli-> 

2 VM 2 M, 
< 01^|0 + > 



AyVsHC + 20> e^l £ 2 ^, ys^ + 2V^> , (* = 1, 10) 



< r|£|0 + > , , , 

5 = y 2 Wo, s* 2 , -s*\ 0)£\ 

<0 "!j!!l > = 7^(0,^-^,0)^ (* = 4,5), 



2/ yjM 2 M x 

< 01/,'|2 + > 
2/ VM 2 Mi 



(A: = 6,7) 



< l~ljg|2 + > 
2 VM 2 Mi 

< Hi1l2 + > 
2 VM 2 Mi 

< ng > 

2 VM 2 Mi 



2j VM 2 Mi 
-1 ;3ii- v 



= r 2 v( r y e 33 + 2 r 2 {^ + *g>] , e 2 ^ 

r£ -[^ + 2 r V { ^-^)] 
= 2 r V 2 ( e ! 3 e f- £ f e * 1 ){^-2 r 2 ^}, 

= 2 r V( e ! 3 4 2 - e ?V 1 ){^ 1 ) -2 r 2 y 2 ^}, 
= 2 r 4 yj (e f e f- e f £ f)^ + v 2 



■(8) 
,4 1 ' 



(£[ &2 + Sj e 2 



.33 .2^(9) 



-/y 



(coe l e 2 + s l s 2 )£ A3 +4y V |£j e 2 t; M - £j e 2 £ 



.(9)1 
'A5 I 



Hill! 



2j VM 2 Mi 

41 



2/ a/M 2 Mi 



(Al) 
(A2) 
(A3) 
(A4) 



(A5) 



(A6) 



(A7) 



k 2 {(« 2 4 3 - e -te 2 )C + V 2 (. 2 4 3 + ^4 2 C} 

-2y^[s\s^-s\sf^), (A8) 

where the repeated roman indices should be understood as contraction with respect to the spatial componets perpendicular to 
the momentum, that is, s\sV = 2/=i,2 £ \ s 2 ■ The lower indices 1 and 2 of M and the polarization vectors/tensors stand for the 
initial and the final one, respectively. 



Appendix B: The angular-spin part of the wave functions 

Polarization vectors and tensors of mesons satisfy the following orthonormal conditions and completeness conditions, 

(g(p> . g(0) = S„ (p,q= 1,2,3), (Bl) 



p 

Tj s T e t= 6 ™ (P.9 = 1.2,3,4,5) 
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(B2) 

(B3) 
(B4) 



The angular-spin part of the wave functions of various spins take the expression with the polarization vectors or tensors as 
listed in Table [VII 

TABLE VI: the angular-spin part of wave functions of low lying states 



2S+, Lj 


k 


yip) 


'So 


-1 


1 


3 S. 


-1 


& ■ 




1 


it ■ & 




1 


(it- <?)((? -gV>) 




-2 


i {3n •£ ( '"-(;t •<?)((?•£ 


3 Pi 


-2 




3 D, 


2 





They satisfy the orthonormal conditions, 



4tt2 L/ y 1 p 



(B5) 



where p and g represent both the polarization state and the spin state. 

Appendix C: The decay widths in terms of the radial wave functions 

The transition amplitudes are expressed by the radial wave functions and the spherical Bessel function, f(x) = sinx/x, as 
follows. 

(1) l"( 3 Si)^ ()-(%) 



€1 = y C + 2 = - j drf'( q r)( U2Vl -v 2Ul ) + 0(V 3 ) 



f(qr)u 2 ui - {f(qr) + 2f"(qr)} v 2 v x 



+ 0(V 2 ) 



AD 
»A1 



dr 



f(qr)u 2 u\ + f'\qr)v 2 v\ - Vf'(qr)(u 2 v\ + v 2 U\) 



+ 0(V 2 ) 



r = 



(gQ 2 k R M 2 



2 \ U ^A\~ k ^Al 



(2) + ( 3 /» ) -> O-C'So) 



$A1 = "J drf(qr)(u 2Vl - v 2 u x ) + 0(V 2 ) 

1? = V ^A2 = ~ J drf'(qr)( U2Ul +v 2Vl ) + 0(V 3 ) 



_ (gQ 2 M 2 k R , (2) (2) 
r " 2^M, 1^!-^ 



(CI) 
(C2) 
(C3) 
(C4) 

(C5) 
(C6) 

(C7) 
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- ! d {- f 



f{qr){u z m - v 2 v0 - Vf(qr)u 2 vi + Vf"(qr)v 2Ul + 0(V 3 ) 



r = 



(4) 1+C^i") -» O-C'So) 



V^ 4) = J <frj-/'( 9 r)("2«i - v 2 vj) - V/fer)«2Vi + Vf"(qr)v 2 m J + C(V 3 ) 



r = o 



(5) lT'W^OtSo) 



= J rfr|/'( ? r)( M2Ml - v 2 v0 - | {/( g r) + 3f"{qr)} « 2 v, 



+|{3/(^)+/"( 9 r)}v 2Ml 



0(V 3 ) 



r = o 



(6) lT^i'^ltSi) 



A6) 

Mi 



= J drf(qr)(u 2 vi-v 2 u l ) + 0(V 2 ) 
nf = y ^A2 = j drf(qr){u 2 m + v 2 v0 + 0(V 3 ) 



Vf® = V^ 6 4 ' = J rfrj-/'(< ? r)( M2 K 1 - 

r = ^^L^-Mf 



(7) iT'pn^i-esi) 



^j=f dr \f(ar) + 3f"(qr)} {u 2Vl - v 2 u,) + 0(V 



-f{qr)u 2 m + Wiqr) + 3f"(qr)} v 2 v, 



+ 0(V 3 ) 



V2/ 



f'{qr)u 2 m + [2f'(qr) + 3f"'{qr)} v 2 vi 



~ {/(r) + 3/"(?r)} u 2Vl - I {3/( 9 r) + 5/" (<?/-)} v 2Ml 



V2/ 



<r/r 



-2f'(qr)u 2 ui - {f'(qr) + 3f"'(qr)} v 2 v\ 



+V {f{qr) + 3f"(qr)} u 2 v x + 2Vf"(qr)v 2U] 



r = 



(gQ 2 k R M 2 1 
27TMJ/2 3[ 



+ <9(V j ) 

2 



- k B r,^f x 2 + [u B r,«l - k B T%} 



T t \ 0}B ^Al~ k B r lA 2 



(C8) 
(C9) 

(CIO) 
(Cll) 



(C12) 
(C13) 

(C14) 
(C15) 
(C16) 

(C17) 

(C18) 
(C19) 
(C20) 
(C21) 

(C22) 
(C23) 
(C24) 
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(8) 2 + ( s P 2 ) -» O-CSo) 



(9) 2 + ( j P 2 )^ 



(8) 



(8) 
%2 



<9) 
A3 



(9) 
%3 



(10) l-^DO^O-CSo) 



^J3 jdr 



V^> + 2^ + 2^} 

J<ir If for) + 3/' V)} (ifcvi - v 2 mi) + 0(V 4 ) 

£(8) , 9v 2 M 6) _ £ (6, 1 
Ml iv IM2 'A3 J 

-f{qr)u 2 m + (2/'for) + 3/"'for)} v 2 n 
-f'(qr)u 2 u x - {/'for) + 2/"'for)} v 2 v x 

+ 0(V 3 ) 



0(V 3 ) 



V 



r = 



(9) 



(9) 
%2 



+ 2 {fiqr) + 3f"(qr)} (u 2 Vi + v 2 m) 



(80 k R M 2 j ( 8) , (8) 



J C /r{/for) + 3/'V)}( M2 v 1 -v 2 M 1 ) + 0(y 



/'for)«2«i - {2 f{qr) + 3/"'for)} v 2 v l 
-{f(,qr)+f"(qr)}v 2 v i 



+ 0(V 3 ) 



+V{f(qr) + f"{qr)}v 2Ul 



+ 0{V S ) 



= 2V 



-£} + V^ + 2£ 



V3 j dr 



-f(qr)u 2 ui - [f'{qr) + 2f"'(qr)} v 2 v\ 



V 



+ 2 + 3 /"(^)) ("2Vl + V 2 Mi) 



(gQ k R M 2 \ (9) (9) 

— TTi^Al 



+ 0(y J ) 



C = V(^ 1 1 0) +2^ A 1 2 0) ) = -^ J Jr/'forX^vi-vzwO + ^V 3 ) 



(10) = -(10) . 7V 2AW 
'lA2 ~ ?A1 " l " z,v M3 



V2/ 



_ ( g ^) 2 ^M 2 r (10) (10) 

r - (mMj* t bT1m ~ kBT,A2 



\f{qr) + 3/" for)} m 2 mi + {/for) - /"for)} v 2 Vi 

2 



+ 0(V 2 ) 



(11) r( 3 D!) -» r( 3 5!) 



n ( M = V(fla + 2 = / drf(qr)(u 2Vl - v 2 m) + 0(V 3 



(C25) 

(C26) 

(C27) 
(C28) 



(C29) 



(C30) 



(C31) 



(C32) 
(C33) 

(C34) 

(C35) 
(C36) 

(C37) 
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(g0 2 k R M 2 [ a,) , (ii)V 



(C38) 



(C39) 



As defined in Section II, q = -2mg V. V is the velocity of the heavy-light meson in the Breit frame. u>b and kg denote the energy 
and momentum of the emitted chiral particle in the Breit frame, respectively. which appears in the final expression of F, is i 
times the coefficient in front of the chiral field in the matrix Eq. 
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